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ABSTRACT 


This report describes improvements and extensions to the RESIST computer 
program, developed previously, for determining the normalized modal 
stress response of Shuttle insulation tiles . The new version of RESIST can 
accommodate primary structure panels with closed-cell stringers , in addition 
to the previously developed capability for treating open-cell stringers. In 
addition, the present version of RESIST numerically solves vibration problems 
several times faster than its predecessor. 

A new digital computer program, titled ARREST (Acoustic Response of 
REusable Shuttle Tiles) is also described. Starting with modal information 
contained on output tapes from RESIST computer runs, ARREST determines 
RMS stresses, deflections and accelerations of Shuttle panels with reusable 
surface insulation tiles o Both programs are applicable to stringer stiffened 
structural panels with or without reusable surface insulation tiles o 
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Section 1 
IWTRODUCTIOW 

This report describes extensions to vibration and stress analysis work 
performed on the Space Shuttle thermal protection system (Reference i)o These 
extensions axe concerned with improving the struetureJ. idealization and 
computational efficiency of the modal response programs developing an 
associated acoustic response analysis capabilityo 

In the initial work, an efficient iterative procedure for the vibration 
and modal stress analysis of Reusable Surface Insulation (RSI) of multi-tiled 
Space Shuttle panels was presented. The method, which is quite general, is 
rapidly convergent and highly useful for this application. A user-oriented 
computer program based upon this procedure and titled RESIST (REusable Surface 
Insulation STresses) was prepared for the analysis of soft-bonded Shuttle 
insulation tiles attached to stringer-stif ened panels. RESIST, which uses 
finite element methods, obtains three-dimeni. ional tile stresses in the soft 
isolator and RSI materials. In addition, the program can treat non-uniform 
temperature-dependent material properties in the tiles. Two-dimensional 
stresses are obtained in the tile coating and the stringer-stiffened primary 
structure plate. A special feature of the program is that all the usual 
detailed finite element grid data is generated internally from a minimum of 
input datao The program can accommodate tile idealizations with up to 850 nodes 
(2550 degrees-of-freedom) per tile and primary structure idealization with 
3000 nodes . The Vibration capability of this program was made feasible through 
development of a new rapid eigenvalue program named ALARM (Automatic LArge 
Reduction of Matrices to tridiagonal form) which can handle primary structure 
idealizations with up to 15,000 degrees-of-freedom. 

Because of the complex nature of the finite element idealizations and the 
large size of the associated matrices. Computer running time was an important 
factor associated with operation of the RESIST program. Therefore, attention 
was given in the present study to identifying time-consuming computational 
areas and developing improved programming and numerical algorithms for certain 
portions of RESIST. These improvements, which are discussed in Section 2.1, 
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resulted in efficiencies which have yielded computer cost-savings of over 75/5 
for realistic tile problems. 

Another basic improvement in the original version of RESIST was the 
development of an improved stiffened panel idealization. The version of 
RESIST described in Reference 1 could only accommodate compact, widely-spaced 
stringers, or open-section stringers attached along a single rivet (or weld) 
line. The present work, as described in Section 2.2, extends the RESIST 
modeling capability to the treatment of orthotropic plate idealizations , or 
closed-cell stringers attached to flat panels along double rivet lines. 

Thus, RESIST has been modified to treat configurations which are more 
representative of actual Shuttle structiore. 

The vibration analysis of RSI tiled panels by the RESIST program yields 
normalized deflections and tile stresses for each mode considered. For dynamic 
analyses, the results for each natural freq.aency may be multiplied by suitable 
modal participation factors and the results superimposed to yield useful 
engineering results. A computer program for obtaining the acoustic response of 
(spatially) uniformly loaded panels has been developed. The new program titled 
ARREST (Acoustic Response of REusable guttle Tiles) is described in Section 3. 

Typical results for the combined usage of ARREST and the new version of 
RESIST are presented in Section E. Preliminary resiilts reveal that the RSI 
normal thickness stresses, near the isolator bond line, may be critical for 
portions of the orbiter which are subject to high acoustic launch loads. 

Users manuals for these programs are presented in Appendices A and B. 

Reference 2 is the Programmer’s Manual for RESIST. Appendix C contains a 
programmer's manual for the newly developed ARREST program. The RESIST and 
ARREST computer programs have been developed on digital magnetic tape for both 
the GDC 6600- and the IBM 3T0-series machines. 
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Section 2 

IMPROVEMENTS TO THE RESIST PROGRAM 

A nurater of changes to the RESIST (RSusahle Siarface Insulation OTresses) 
program (References l and 2) have heen made to improve its usefulness. The 
nature of these program changes have heen to extend its applicability, reduce 
the program storage req.uireraents , speed up the machine running time, and to 
provide additional information which is required for subsequent dynamic 
analyses (described in Section 3). 

2.1 COMPUTATIONAL EFFICIENCIES 

The reduced program storage requirements and the development of streami- 
lined numerical algorithms have improved the program efficiency significantly. 

This is highlighted by a reduction in machine running time of a large sample 
problem, presented in Reference 1, from 68 to l6 CPU minutes on the IBM 370/168 
computer. These program improvements will now be detailed. 

2.1.1 Modification of The Large Eigenvalue Routine (ALARM) 

The ALARM (Reference 1) routine computes the initial mode shapes and fre- 
quencies of the tiled plate. A modification to ALARM was effected to speed up 
this time consuming portion of the RESIST computer program. In the original 
version of ALARM, the number of eigenvalue reduction vectors (Reference 3) gen- 
erated, {Vj}, 1, 2. ..., ra, was approximately twice the number of desired 

modes i= 1, 2, ..., q. Since the generation of each Vj represents one 

of the most computationally costly portions of the ALARM algorithm, the sub- 
routine logic was changed to initially compute only approximately 50^ more 
reduction vectors than desired eigenvectors (m « 3/2 q) . The first q eigen- 
values associated with the m reduction vectors are now computed along with 
their associated Newman-Pipano error bounds. Should any of the error bounds 
in frequency squared be greater than 2^*, an additional vector is computed and 
the first q eigenvalues and error bounds, using m + 1 reduction vectors this time, 
are checked again. This process is repeated, adding one new vector each step 

*An error bound in frequency-squared of 25? corresponds to one of approximately 
1/? on frequency. 
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of the way, until the first q. desired eigenvalues are accurate to within 25?. 

It has been fovmd that this procedure is e^ctremely efficient. Repeated solu- 
tion of the reduced m x m eigenvalue problem, and error bound deterpiination , 
requires a minute portion of the entire computation time , as long as 5 or more 
modes are required and the size of the original eigenvalue problem is large 
(i.e., of the order of 500 or more). 

2.1.2 Monpositive Definite Solution Algorithm 

The efficient solution routine (PODSYM), used to solve the tile and 
primary structure banded large matrix equations, required that the symmetric 
coefficient matrices be positive definite. Since this condition could not be 
guaranteed a priori, Eqs. ( 6 ) and (ll) of Reference 1 were converted to the 
equivalent positive definite forms of Eqs. ( 9 ) and (12), respect- -ely. This 
conversion required excessive storage requirements and resulted in rather high 
computer usage costs. Therefore, PODSYM was modified to accommodate nonpositive 
definite equation systems through decomposition of the symmetric coefficient 
matrix [A] into the form 


[A] = [L] ['DJ [if 


where [L] is a lower triangular matrix with elements 1. . and [’’D.,] is a 
diagonal matrix with plus or minus ones , d , on the diagonal . The original 
algorithm is similar to the one outlined below except that the computation and 
inclusion of d , which is always plus one if [A] is positive definite, is not 
necessary. The elements of [L] are given by 


Z. . 

IX 



i -1 


■-S 


(“ii - 2 %) 

p=i 


0-1 

‘ij - 2 

p=i 


a . . - > d Z, I . 

xj P xp jp 



2-2 



where 


i-i 

d. = sign {a. . - 7 d ) 

1 ® 11 p xp 

p=l 

and a., are the elements of [A], 
ij 

The solution of 
[A] {x} = {b} 

for {x} (given {b} ) then follows from the forward solution of 

[L] {y} = {b} 


for {y},and the backward solution of {x} from 
[L]'^ {x} = ["Lj {y} . 


2.1.3 Iterative Procedm*e for CSianging Matrix Equations 

The iterative solution of Eqs. (6) and (ll) of Beference 1, both of which 
are of the form 


[Z.] = {p (5.) } 


( 1 ) 


where 


[Z.] = [K] - 0,2 [M] 

required the repetitive decomposition of [Z^] into [L. ] [^D.,] [L^] , since 
the frequency, upon which [Z^] depends, changed from one iteration to the 
next. However, since o,^ may be a slowly changing value, an option, which is 
based upon an approximation, and which does not require the resplitting of 

was introduced into EESIST. This approximahion is based upon the follow- 
ing development. 
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(U^ = (x)^ + Aw? 

1 O 1 




' «Iil> " 


and substitute into Eq. (i) to obtain 




( 2 ) 


where the higher order terms Aw^ [M] l^a-ve been dropped from the right- 

hand side of Eq,. (2). Starting with = 6^^, Eq. (2) is solved in an iterative 

manner (this is an "inner" (j) iteration) until converges. Thus, 

[Z^] is only decomposed once, regardless of the number of "outer" (i) 
iterations . 

2.2 PROGRAM EXTENSIONS 

The extended applicability associated with the new version of RESIST 
involves improvement of the idealization for the stiff ened-plate primary 
structure. The new plate idealizations permit the use of either orthotropic 
plate properties or a more detailed stringer model. The new finite element 
capability for the strain isolator includes the option for material orthotropy, 
which is required for the current baseline TPS isolator. 

2.2.1 Orthotropic Strain Isolator 

RESIST has been modified to permit an orthotropic material isolator 
between the RSI and primary structure. This latter modification was necessary 
because of the current baseline usage of Nomex felt (which has different 
properties in the transverse and planar direct! 'vs) as the isolator material. 

The finite element modeling of the isolator layer is still based upon 

the three-dimensional orthotropic hexahedron element described in Reference 1. 


2-h 


2.2.2 Stiffened Plate Idealization 


The stiffened-plate primary structure finite-e].ement idealization, 
employed in the original version of RESIST, consisted of offset beam elements 
attached to flat plate elements (Figure 2-1). While it is felt that this 
idealization is adequate for compact widely-spaced long stringers, it is 
inappropriate for a "typical” Shu’ tie panel, such as that shown in Figure 2-2, 

T: oi.teiK- the capability of RESIST for application to such structural con- 
I'ifuratior.s, two new idealization options have been developed and incorporated 
into tnis program. 

The first development consists of a new stringer idealization co permit a 
more accurate modeling for closed-section stringers attached along a double 
line of rivets (or welds). The new stringer model consists of a pair of plate 
elements, which form a "V" (Figure 2-3) > below which an offset beam element is 
rigidly attached. The rationale used in selecting the "V" shape was that it 
permitted the modeling of a shear-flow path to resist local plate twisting (the 
previous offset beam idealization could not furnish this type of twisting re- 
straint) while only adding a single row of element-nodes per stiffener. 

The offset beam element placed below the "V" is used to compensate for 
the fact that the original stringer possesses more cross sectional beam prop- 
erties than can be accommodated by the "V" elements alone . Thus , when an actual 
closed stringer to be idealized has a shape that is different from the one 
used in thi- finite element model, the program assigns properties to the offset 
beam below the "V" such that its -.<.-ctlon properties compensate for any mismatches. 
A count of the actual n^iraber of stringer sect.ion properties yields six geometric 
values (a, 1,1, , and e ) in addition to the spacing variables P, a and b 

(Figure 2-2), whereas the new idealloa:: ion permits seven independent properties 
in addition F, a ana b. The arlditloual variable permitted by the present 
finite element idealization, T, may thus be used to help impjrove correlation 
with tent lata. 

The second sti.’uct-ural idealization option incorporated within RESIST 
permits the i.-ombination of stringer properties into the flat plate elements. 

This in affected in a uniformly distributed manner through specification of 
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Figur« 2-1 Stiffened-Piate Finite Etement Idealization Employed in Original Venion of RESIST Program 
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Figure 2-2 Typical Shuttle Panel 
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Figure 2-3 Idealization of New Stringer in RESIST Program 




appropriate orthotropic plate properties. Additional plate input required to 
describe the orthotropy includes the elastic constants Ai 2 s ^22 A 23 

where 



The general form of the above stress-strain law allows tb^ use of various ortho- 
tropic plate theories depending upon the particular specification of the Aij- 
The equivalent values of these elastic constants which were suggested by 
Timoshenko (Reference 5) and Fung (sej Figure 2-h and Reference 6) are compared 
with their isotropic counterparts in Table 2-1. The form of the governing ortho- 
tropic plate equations for these four theories is presented in Table 2-2. 

2 . 2.3 Additional Modal Information 

The new version of RESIST has been extended to output necessary input 
information for the ARREST program described in Section 3. The new version of 
RESIST computes tae modal area terms required for Eqs. (l3) and (l4) of Section 3- 
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Table 2-1 

Comparison of Orthotropic Plate Variables 
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Table 2-2 

Comparison of Orthotropic Plate Equations 














These quantities are related to M. through the relationship 



(i) 

PS 



RESIST has also heen modified to compute norraEillzed modal stringer stresses 
along with the RSI and plate stresses. This data is outpjut on magnetic tape, 
together with modal deflection data required for the acoustic response 
program. 
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Section 3 

SHUTTLE TILE RESPONSE TO RANDOM PRESSURE FLUCTUATIONS 

A computer program for calculation of tlie response of shuttle panels with 
surface insulation subjected to uniform, spatially~Gorrelated acoustic excita- 
tion has been developed, This program, titled ARREST (Acoustic Response of 
REusable Shuttle ^iles) was designed to interface with RESISTo For a dynamics 
problem, the output of RESIST (which consists of modal frequencies and normal- 
ized stresses of stiffened panels with RSI tiles) serves as the input to ARREST, 

The output of ARREST is approximate in that the analysis is based upon a number 
of simplifying assumptions. As will be seen, these assumptions are the usual 
ones associated with the use of finite element analyses, modal methods, and a 
simplified characterization of the acoustic forcing function. Perhaps the most 
limiting approximation included within ARREST is that the normal pressure field 
is uniformly correlated over the entire structural panel. This is a conservative 
assumption that is justifiable for excitation where the pressure field correlation- 
length is of the order of the panel size. When the excitation is caused hy jet 
noise, such as dtiring Shuttle lif"ioff , the correlation length is of the order 
of the engine diameter (Reference T). The theory upon which ARREST is based is 
summarized below, 

3.1 MODAL DYNAMICS OF RSI TILES AND PRIMARY STRUCTURE 

The finite element equations governing the dynamic response of RSI tiled 
panels are of the general form 
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v/ltere {SJ- ai-e xhe I'Odai degi-ees-of- freedom; [Ml, [C] and [k] are the mass, 
viamping and stiffness waT.rices , respectively; and {P} is the dyuajnic loadJc;', 
vector. The snbscrip^n "PS" and "T" denote the primary structvire and tile 
components, j'especti vely. The significance of the bar over the tile quantiti.er; 
Is that they refer t-i the accLumolation of all tiles rather than a single tile, 


It shooj.d be noted rhat the primary structure and tile matrices have beer, 
oast in imccupj.ed fom. The reason for this is that the iterative solution 
procedui-e employed to obtain the matrices and modal vectors of the coupled 
syotem vas to t,-;v.-ai,- ertoh component separately. Their initial coupling is crci- 
itainei vithan r.be fncT. xhat the tile/primary structui-e interfa.ce degrees-nf- 
freedom, a subset of itp, 1 an'3. {6^}, are kinematically dependent, as are tiic 
vectors which contain their mutual reactions. Therefore, 

elements of which correspond to {P^} have, by Newton's Tirird Lavr, the 

same magnitudes but ox>P'^aite signs. 

The classical modal approximation is now made; 
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where the number of orthogonal modes, y 
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and generalized coordinates 


, are less tluirj the number of nodal degrees-of-freedom. The modal vsctoi-s 
satisfy the maxrlx equations 
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where (P.) are the tile boundary interface loads acting upon the primary 

iJ 

structure and {-P } are their reactions, which act upon the corresponding tile 
J 

nodes. Thus, the degree-of~ freedom magnitudes at these common nodes will be 

(i) “(i) 

identical in the vector components r and {(})„ /. 

The modal normalization condition selected Is 


^PS ^PS 


(5) 


and the orthogonality of the modes yields 


■^^PS I 



’“ps i 0 ’ 


° 1 \ 

Ti’*'! 



L^t j 


1 b-j 


where M. = [M^g] [^] 


( 6 ) 


is the modal mass and 5^ ^ is the Kronecker delta with properties 


I = 1 for i = j 

Y 

1=0 for i ^ j 

Substitution of Eq.. (2) into Eq,. (l) , premultiplication by 

u<i> ; and application of Eqs. (3)» (^) and (6) yields the scaler 

PS I -L 

equations 

H. (q. + 2 ? w 4 + q. ) = Q. 

11 111 xl 1 


(?) 


where the generalized forces, Q^, are given by 


Qi= iPps>+ Ui^^J {P^> 


:(i) 


( 8 ) 
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_The unlrnown damping matrix is assxmied to yield uncoupled modal damping ratios, 
(for which values are either assiimed or are measured in test), where 

2 C. M.w. 6., = [C ] 

11 ij ^ps ps ^ps 

+ U^J [Cji (9) 

and u^e has "been made of the relationship 

* Li<">J tKjl a'^h (10) 


3.2 ACOUSTIC RESPONSE TO UNIFORMLY CORRELATED PRESSURES 


If the forcing function is a stationary, random, acoustic pressure and 
the system damping is very small, the mean square response, in the m^^ 
degree-of- freedom, 6^, may be approximated as (References J, 8) 



m 


1 11 


S. . 
11 


( 11 ) 


S^^ is related to the pressure power spectral density, [S (u^)] throiigh. 
the relationship 


Sii == La [Sp(oi.)l {A (12) 

where [S^(ai^)] (measured in psi^/radian) may vary spatially over the structure, 

and {a ^ } are the modal degrees-of-freeciom multiplied by suitable terms 
which account for the local normal surface area over which the pressure 


excitation acts . 


If the pressure distribution is assumed to be "uniformly correlated" 
(i.e., the time averaged mean square random pressure does not vary spatially 
over the structure) then|Sp(a)^ )J is a matrix with equal elements S(u^) and 
Eq. (12) simplifies to 


S. . 


XI 


S(uj^) 



(13) 


where (j)^ are the tile degrees-of~freedom, in the i mode, which are 

l,Zj 

normal to the pressure loading and are the corresponding areas over which 
the pressure loading acts . 


.th 


A similar approximation for the stress component of a given tile at the 


j node, under uniformly correlated random pressvire loading, yields 


,(i)^ 


=2 TT 


i 11 


S(w^) 




ilh) 


for the mean square tile stress. It should be noted that the coefficient for 
Eqs. (ll) and (ll+) changes from tt/ 4 to 1/8 if S is given in psi^/Kz, i-eo. 


f s<“i> =is(f.) 


"til 

where f^ is the i modal natural frequency in Hz. 
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Section 4 
NUMERICAL RESULTS 

Finite element results were obtained for an l8x54-inch stringer-stiffened 
panel, with and without RSI tiles. Three different structural models were 
employed to idealize the panel without tiles. Based on modal results from 
RESIST, the most appropriate of these was then used to model the panel with 
tiles. Modal and acoustic response results were then obtained for an assumed 
loading spectrum associated with an aft section of the Shuttle orhiter fuselage 
at launch, 

4.1 UNTILED PANELS 

The three stiffened-panel options included in the new version of RESIST 
(Section 2.2.2) were used to analyze the Shuttle test panel of Figure 2-2. 

A comparison of the detailed stringer idealizations, with the actual cross 
section, is shown in Figure 4-1. 

Partial frequency and mode shape results for the three models are 
presented in Table 4-1 and Figures 4^2 and 4-3. These results indicate that 
the structural stringer model selected is not critical for the "beam" type 
modes (m=0), but that the cross-stringer modes (m > l) are sensitive to the 
stringer idealization, with the orthotropic (Figure 4-lb) and closed-cell 
stringer (Figure 4-ld) options agreeing more closely with one another, than 
the single attachment line stringers (Figure 4-lc). This is not surprising 
as the m > o modes require the transmission of torsional load through differ- 
ential shear, which the singly-attached stringer cannot model properly, whereas 
the m=0 modes primarily involve stringer bending effects. 

4.2 TILED PANEL RESULTS 

The RSI-tiled panel of Figure 4-4 was analyzed by the RESIST program using 
the stringer idealization of Figure 4-ld. There are l8 6x6-inch tiles between 
supports, each of which is 2.6 inches thick. There is a 0.01-ineh thick 
coating over the top of each tile and the isolator thickness is 0.16 inch. 
Nine-inch extensions, over opposite edges of the supports, were assumed to act 
as rigid bodies with rotatory inertia properties comparable to the interior 
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Figure 4-1 Panel Stringer Idealizations Used for the Shuttle Panel (Fig. 2-2) with the RESIST Program 


Table 4—1 Natural Frequencies (Hz) fur the Shuttle Panel (Figure 2-2) Computed by RESIST 

for Three Different Stringer Options 


m* 

n** 

Stringer Attached 
to Plate Along 
Single Rivet Line 

Closed-Cell Stringers 
Attached Along 
Double Rivet Lines 

t = .008 in. 

Orthotropic 

Plate 

Option 

Dx = 220,000 lb-in. 
Dy 31 lb-in. 

Dxy ' 9850 lb-in. 

0 

1 

106 

110 

107 

1 

1 

108 

133 

134 

2 

1 

126 

174 

184 

0 

2 

286 

285 

271 


*m = number of nodes in cross-stringer direction 

**n ~ number of % sine waves between spanwise supports 
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Figure 4-4 Idealized Shuttle Panel (36 x “IS Inches) with 18 RSI Tiles 
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(36-inch) section of tiled panel. The finite element idealization of a typical 
tile is shown in Figure k-5 . The resulting natural frequencies and modal data 
are given in Table k~2. 

The random sound pressure level curve used to simulate the panel's 
excitation (Figure i*-6) is thought to be representative of a structural region 
on the Shuttle orbiter base at liftoff. Detailed results are given for 3% 
modal damping. However, some results for 1% of critical damping are also 
presented in Table 4-3 because of the critical nature of the Shuttle TPS and 
the uncertainty associated with the damping parameter. 


Since the AEREST program assumes uniform spatially-correlated random 
pressiure, and the structure contains two planes of symmetry, only certain 
modes will contribute to the acoustic response. The contributing modes 
correspond to those for which the number of spanwise half-waves (n) is odd, 
and the niamber of chordwise-noles (m) is either 0 or even. This effect is 

reflected by the modal termsf / A . of Table 4-2 (see Eq. (ll) and (l4) 

j f 

of Section 3). As may be seen, this term decays rapidly beyond the funda- 
mental, even for the symmetric modes. Thus, for most practical purposes, 
only the first couple of modes need be considered. 

Uorraal-to-the-plate RMS acceleration and stress levels for this excitation 
are presented in Figures 4-7 and 4-8, As determined by an automatic search 
procedure of the results, which is incorporated within ARREST, the o tile 
stresses are highest. In addition, the lowest RSI Material allowables are also 
lowest for these stress components. Hence, only the a stresses are plotted 
in Figures 4-9 through 4-13. Smooth curves are not used for the tile RMS stress- 
plots as the plan-form idealization is composed of a finite element grid that 
is comparatively coarse (4x5), hut adequate, for the present purposes. The RSI 
stresses are plotted for Z = .485 iu. since these are the specific points within 
the highest stressed elements for which the stress components were actually 
calculated. As shown in Table 4-3, the peak stress levels become critical for a 
scatter-factor* of 3. This is especially true when the modal damping decreases 
to 1%, 

*A scatter factor of 3, applied to the RMS value, simply means that the 
probability that a given stress component will be below this value at 
any instant is approximately 99 o 5^. 
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Figure 4-5 Finite Element Ideall^tion of Typical RSI Tiles 
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Table 4-2 

Computed Modal Data for the Shuttle Panel 
with 2.6-inch RSI Tiles (Figure 4-4) 


Mode 

Number 

i 

MODE 

m* 

SHAPE 

n** 

Frequency 

Hz 

(i) 2 

J ' * 

(Refrirence Eqs. 13 and 14. 
Section 3) 

Sound 

Pressure 

Level 

S(f.) 

(psi^/Hz) 

1 

0 

1 

72 

627.5 


2 

1 

1 

81 

0 


3 

2 

1 

92 

1.186 


4 

3 

1 

102 

0 

HI 

*m = number of nodes in cross-stringer direction 
**n = number of % waves between spanwise supports 
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Figure 4-12 Critical Acoustic Stresses Aaoss Tiles Near Supported-Edge of Panel 








table 4-3 

Critical Tile Stresses for Acoustic Launch Loads 


Tile No. 

(See Figure 4-2) 

r-inite 

Element No. 
(See Figure 4-S) 

azi(psi)* 

3% Modal Damping 

1% Modal Damping 

RMS 

3a** 

RMS 

3a** 

7 

33 

2.317 

6.951 

• 4.013 

12.039 

7 

28 

2.203 

6.609 

3.816 

11.448 

9 

33 

2.150 

6.450 

3.724 

11.181 

9 

28 

2.094 

6.282 

3.627 

10.880 


*The allowable Ozz for the RSI at room temperature is approximately 8 psi 
**The 3o values correspond to a scatter factor of 3 (i.e., approximately 99.5% probability that the 
instantaneous stresses will be below this number) 
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Section 5 

CONCLUSIONS AND RECOMMENDATIONS 

A user-oriented program, RESIST, has been developed to compute the static, 
thermal, and normalized-flaodal stresses of arbitrarily heated RSI tiles affixed 
to typical Shuttle structural panels. The current phase of this effort reports 
upon improvements to the modal portion of RESIST and the development of a new 
program, ARREST. ARREST, which computes the response of the tiles and panels to 
uniformly-correlated random acoustic pressures, was designed to efficiently 
interface with the RESIST program. 

The present version of the RESIST program has been streamlined and runs 
about four times as fast as its predecessor. In addition, it has been extended 
to accommodate orthotropic plates or isotropic plates with tmiformly spaced, 
arbitrarily shaped stringers which are either singly or doubly attached to the 
plate. It is believed that the orthotropic and doubly attached stringer 
options better represent the stringer twist-stiffening, which is required to 
accurately predict the plate cross-stringer modes (m > o). 

Results from the ARREST computer program for a simulated Shuttle launch 
environment indicate potentially critical normal thickness stresses in the RSI. 
The other RSI stress components are generally lower and possess strength 
allowables which are 3 or it times as high as well. 
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Section 6 
NOMENCLATURE 


[*] 

[c] 

bi 

b] 

b] 

[m] 

pp] 

pi] 

p] 

IM 

fp| 

PI 

PjI 

1*1 

M 

Pil 


lP^>t 




E 


Symmetric matrix 
Damping matrix 

Diagonal raatri i with plus or minus ones on diagonal 

Stiffness matrix 

Lower triangular matrix 

Mass matrix 

Pressure powei’ spectral density matrix 
Impedance matrix associated with 
Modal matrix 

Right hand side of matrix equation [a] |x| = {hj 

Forcing function vector 

Generalized modal response vector 
*fe V| 

J eigenvalue reduction vector 
Solution vector of equation [a] {x| = {lil 
Forward solution of matrix equation [l] {y} = Ibj 
Nodal displacements associated with i iteration 

i modal vector 
Stringer cross-sectional area 
Orthotropic plate stiffness properties 
Area normal to z direction at j node 
Plate hending stiffness 
Plate modulus of elasticity 



E’, E', E" 
x’ y’ 




S (to.) 
a, b, P 


d 

P 

e 

t 

h 

m 

q. 

t 


w 


X, y, z 
xV, y’ , z’ 

e 


Elastic constants associated with orthotropic plate theory 
(see Table 2-1) 

Plate shear modulus 

Stringer moment of inertiaL about local principal axes 

Stringer polar moment of inertia about principal axis 
th 

Modal mass for i mode 

Primary structure portion of i^^ modal mass 
Generalized modal force 

"fell 

Acoustic sound pressure level at i^ frequency 
Stringer spacing parameters (see Figure 2-2) 

Element of [a] matrix 

Plus or minus one 

Stringer eccentricity 

Frequency in Hz 

Plate thicloiess 

Humber of reduction vectors 

Number of modal vectors 

Thickness parameter associated with stringer idealization 
(see Figure 2-3) 

Normal plate displacement 

Cartesian plate axes 

Stringer principal axes 

Angle between stringer principal axes and plate global axes 


6 . . 

10 




y’ w 


Kronecker delta 
Plate membrane strains 

Critical damping ratio associated with i * mode 


Plate stress components 

RSI direct stress component in Tihicicness direction 

Normal (z) deflection of j tile node associated with 
mode shape 

*th. 

Natural frequency of i mode in radians per second 
Stringer pitch parameter (see Figure 2-k) 

Eleraert of [l] matrix 
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APPENDIX A 
USER’S MANUAL FOR 
RE*S*I*ST 

(STATIC AND DYNAI-IIC REUSABLE SURFACE DSSULATION STRESS PROGRAM) 


AUGUST 1975 VERSION 


PRECEDING PAGE BLANK NOI nUOIl 


A. IKTEODUCTION 


This Appendix describes the use of a finite element based structural com- 
puter program for deteitnining the static response and natural vibrations of 
TPS protected shuttle panels. The program is titled "EESIST” for static and 
dynamic ^usable Surface Insulation Stresses. The logic flow for RESIST is 
presented in Figure A-1. 


The basis for the method is that the TPS is nonstructural but its stress 
levels, which are critical, must be computed. Thus, it becomes possible to 
neglect the stiffness of the TPS initially, but not its mass in the vibration, 
to deteimine approximate primary structure deflections. 

An iterative procedure is then perfoimed whei;e, for each step, the primar'* 
structure deflections are imposed individually upon each tile at the tile/ 
primary- structure interface, and the tile deflections and interface boundary 
loads are obtained. For the vibration option, the frequency is updated by com- 
puting a Rayleigh quotient, using the latest non-rigid tile displacements in 
addition to the corresponding primary structure displacements. The individual 
tile boundary loads obtained are then assembled and their reactions are applied 
to the primary structure. New primary-structure deflections are obtained and 
compared to the previous set. This process is repeated until convergence is 
obtained. 

B. PROGRAM LIMITATIONS 

The usual assumptions for programs based upon the linear elastic finite 
element method are applicable to RESIST. However, - to facilitate the prepara- 
tion of program input, a number of simplifications regarding the configura- 
tion and loadings have been made. Thus, the generation of a voluminous quan- 
tity of finite element input data has been greatly reduced by inclusion of 
a series of data preprocessing subroutines within RESIST. The restrictions 
upon which these subroutines are based follows: 

1, Boundary conditions and edge loadings are assumed uniform along 
the four rectangular plate edges defined by x = 0, t and y - 0, 

A 

V 

2. The primary structure plate temperature and properties are all 
uniform. 


A-1 


ORIGINAE PAGE IS 
OP POOR QUALITY 


m 


INPUT 


IX 


/ 

MESH 


FORM K 

1 * 

GENERATION 


8. SET i = 1 

< VIII 

X 

1 

1 


SOLVE FOR 
PRIMARY 
STRUCTURE 
DISPLACEMENTS 


FORM loading 

OFPRIMARY 

STRUCTURE** 


STATICS 
ITERAtlON I 

*- 


COMPUTE tile 
THERMAL & 
external MECH 
loads 



COMPUTE 

4 VB^ 

FORM 

MASS 

QUOtlENt 


MATRICES 


VI (HEX) 

V lALARMI 1 

f 

OBTAIN TILE 


COMPUTE APPROX 

DEFLECTIONS, 


FREQ & MODE 

MODAL MASS & 


SHAPES 


bound LOADS' 


ZEROTH ITERATION 



APPLY FOAM/PRIMARY-STRUCTUBE BOUNDARY DISPLACEMENT TO 
EACH TILE SPECI FIEO. ON E AT A TIME, AND COMPUTE ALL NON 
boundary Tl LE OISPLACEMENTS AND STEADY STATE BOUNDARY 
LOADS. 

* • FOR VB PROBLEMS USE LOADS FROM VI AND FREQUENCY FROM 
VII. FOR St Problems WITH i = i, us| loads from module x. 
HOWEVER. IF i p 1 GET StATIG TILE LOADS FROM VI AND BOUN. 
OARY LOADS FROM X. 


Figure A-1 Flow Chart for RSI Stresi Analysii Program "RESIST' 
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3. The stringers are equally spaced with temperatures and properties 
which are all uniform. 

U. All tiles are geometrically Identical as are their temperature 
distributions and uniform pressure loadings. 

5. The boundary conditions must be selected such that the primary 
structure is statically stable. 

The remaining limitations are primarily concerned with the program's 
capacity and should be adhered to by the user. These limitations are as 

follows : 

6. Maximum number of nodes in a tile = B 50 . 

T. Maximum number of finite elements running in any one direction 
in a tile =20. 

8. Maximum number of nodes in primary structure = 3000. 

9 . Maximum manber of primary structure nodes along x or y direction = 

1 , 000 . 

10. Maximum number of degrees of freedom in primary structure = 15,000. 

11. Maxiraimi number of natural mode shapes =50. 

12. Maximum number of stringers = 15 • 

A violation of restrictions 6-12, inclusive, will cause the program to 
stop and an appropriate warning message to appear. 

To insure symmetry of solutions for iianels which are symmetric with re- 
gard to stringer locations about y = I. /2, care should be taken with the input 
data to see that the plate nodes associated with the stringers are symmetric 

about y = L /2. 

y 
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C. INPUT INSTRUCTIONS 

A description of the card input for the IBM 370 and CDC 66OO versions of this 
program is presented in this section. 

In addition to the first two input cards which contain literal data, such as 
special program title and date, in columns 1 through 80, inclusive; there are six 
groups of input cards containing the following information; 

• Group A - Instructions regarding the type of problem being performed, num- 
ber of iteratious desired, and type of output information 

• Group B - Details of the geometric configuration and finite element mesh of 
the primary structure and tiles, (Card B,4 is omitted if there are no 
tiles ) 

• Group G - Defines the primary structure and ESI temperature dependent 
material properties. If there is no TPS, cards C,3 through C,ll are omitted, 

• Group D - Specifies the primary structure boundary conditions 

• Group E - Describes the mechanical loading upon the primary structure as well 
as its temperature. These cards are omitted when the vibration option is 
used 

m Group P - Defines the RSI temperature distribution. These cards are 
omitted if there is no TPS. 
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If there are no tiles then n^ and n^, together with n^g and njjg, are still required since they determine 
the primary structure finite element grid. In analyzing panels without tiles, leave out cards B.4, C.3 
through C.IO and all "F" cards. 
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Fig. A-2 Doubly-Connected Stringer Idealiiation - Global Coordinates 
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Figure A-7 Possible Static Mechanical Loadings Upon Panel - Global Coordinates 
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Seven temperature differences to a card until all nodes are accounted for 
Cols. 71-8Q of each card are reserved for user’s card identification. 









D. DESGRIPTION OF OUTPUT 


Output from a typical run of the .-JSIST computer program is 

explained below in outline form. References in parentheses refer to 

pages in this Appendix. 

1. Program title and date indicating latest update of program 
version which was run. 

INPUT IHFORMATIOM 

2. Listing of input cards , the first two of which are the title 
assigned to any given run by the user. 

3. User selected input options are listed (pp. A-5 - A-8) . 

4. Plate, stringer and tile geometry and specification of finite 
element grids for primary structure and tiles (pp. A-9 - A-ll) . 

5. Plate, stringer, strain isolator and arrestor mate riul properties 

(pp* A-1^ - A-17)* Note, if there is no strain arrestor, RSI or isolator 
material properties may be used for the arrestor. If this is done, 
the thickness dimension of the usual isolator or RSI should be 
appropriately reduced to compensate for this addition. 

6. Temperature -dependent RSI material property data used for generating 
curves used internally by program to compute RSI average finite 
element properties (pp. A-l8 - A-19). 

7. Plate and stringer boundary ,c_Qndit ions (pp . A-20 - A-25 ) . 

8. Applied primary structure static mechanical and therma.1 loading 
if not a vibration problem (A-26 - A-^27 ) . 

9* RSI temperatm^e distribution input data. Used for property data (item 6 
above) and thermal loading if a statics problem (pp. A-28 - A-29). 

OUTPUT INFORMTION 

I 

10. showing typical tiles three dimensional finite element ordering, 

by layers. Top, or first layer also corresponds to two-dimensional 
tile coating elements as well. 


A- 30 


11. Map showing ordering of a typical tiles finite element nodes 
by layers. 

12 , Position and temperatures for a typical tile in a local coordinate 
system (reference Figure A-4, A and B). 


13. Global geometry of primary structure nodes and plate nodal degree-of- 

freedom nunibering. D , D and D refer to nodal deflections , and R and 
E are the nodal rotations. Nodes with no degrees-of- freedom are used 

y 

to define the stringer centroids for singly-attached stringers only. 


l4.a. Statics Option: Primary structure nodal deflections by iteration 

number. Nodes with the same x coordinate are grouped together. 

These groups are separated with dashed lines. 

14. b. Vibration Option: Mode numbers, approximate frequencies and 

corresponding modal error bound (which should be less than 2^ 
to be a reliable approximate mode). This is followed by the 
primary structure mode shapes with a similar nodal deflection for- 
mat as for the Statics Option. 

15. If requested by the user, the computed c gnyergenee parameter 

is printed out along with the input quantity it was tested against. 

This is done for each iteration after the first for the Static 
Option. The primary structure degree-of -freedom with the largest 
change from the previous iteration is also identified. 

16. Tile nodal displacements by tile and iteration number. For a 
vibration option, this calculation and the subsequent ones are 
performed only for the user-specified vibration mode. 

17. Modal mass associated with each tile for given mode is printed out 
as DEN. 

18. Three dimensional tile stresses and strains for the bottom two 
layers of elements by element number. These quantities are 
computed at each element's 8 Gauss integration points. Gauss point 
stresses are believed to be more accurate than nodal values and provide 
more detail than simply the element's average stresses. 
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19- Three dimensional element average streases and. strains (by ti3e 
and iteration number). 


20. Two-dimensional element average coating stresses. Coating 
element numbers correspond to three dimensional element numbering 
directly below them. 

21. a. Statics Option: Repeat of items 16-19 for each tile. Repeat 

of item llf.a and 15 for each iteration. 

21. b. Vibrations Option: Computation of Rayleigh Quotient (OMEGA 

SQUARED) if all tiles have been treated. Repeat of items 16-19 
for each iteration. Repeat of items l^.b,15 and Rayleigh 
Quotient until convergence or last iteration is performed. 


22. Plate element stresses and strains for mid and/or top and/or 

■bottom surfaces. This computation is done after each iteration 
if requested by the user. Otherwise, it is computed only after 
convergence or the last iteration is performed. 


23. 


Vibrations Option: Stringer strains and stresses if requested 

in input, the quantity (/ A (printed out as 

*5^ Z ♦ X 5 Z j 


SUM (AREA .* DZ)** 2=), and the primary structure modal mass 
(printed out as P.S. M = ). 


Ar-3? 


13 . SAMPLE PROBLEMS 


Output for a sample statics and vibrations problem, are presented in 
the remaining pages of this Appendix. Only portions of the output are shown. 
However, the pages presertei are representative of the types of information, 
and their respective foi^-ma^s, which the RESIST Program can deliver. 
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APPENDIX B 

USER’S MANUAL FOR 
A*R*RE*S*T 

(Acoustic Response of REusat>le Shuttle Tiles) 


This Appendix describes the input and output of the structural acoustic 
response program ARREST. It is designed for operation with the RESIST computer 
program discussed in Appendix A. The program computes the RMS stresses and 
deflections of stiffened rectangular panels, with or without RSI shuttle tiles, 
subject to uniform correlated (spatially) random acoustic pressures. When 
there ^e TPS tiles, the RSI tile stresses and deflections are edso computed. 

The computation procedure is based upon linear modal methods and assumes well- 
separated lightly damped modes ^ 

The program card input effort is small compared with that for the RESIST 
program, as most of the input is on tapes generated by the RESIST program. 

Each time information for ARREST is generated by RESIST, it is contained on 
tape as units 1 and ,21. Unit 21 contains modal data and unit 1 contains 
structural idealization data. Since the idealization data is the same for 
each mode, only one unit 1 need be input in ARREST. This tape should be 
mounted as unit 9 for an ARREST run. The tapes containing modal data should 
be mounted, in succession, as units 1, 2, 3, T and 8. A minimum of one tape 
with one mode and a maximum of 6 modes may be used in a given ARREST run. If 
information for more than one mode is being combined on unit 1, then unit 1 must 
be mounted with a ring and the added modes are assigned unit numbers 2, 3» etc. 
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^Denoted as (SUM { AREA*DZ ) )**2 on RESIST run output. This g.uantity is printed out at the end of the 
RESIST runs. 


















IN'POT DATA FOB ARREST - Sheet 4 of 4 



**Den.oted as DEN on RESIST run output. This quantity is printed out after eaeh tile displacement. 
***Denoted as P.S.M on RESIST run output. This quantity is printed out at the end of the RESIST runs. 




APPENDIX C 


PROGRAMMERS MANUAL FOR 
A*R»RE»S*T 

(Acoustic Response of REusable Shuttle Tiles) 



• C.l INTRODUCTION 


This Appendix contains the overall flow chart of the program, data set 
allocation and subprogram calling sequence, along with a brief description of 
the individual subprograms and typical subprogram output for the AEffiEST 
program described in Section 4. Two versions of the program exist, one of 
which is written completely in FORTRAN. This version is compatible with 
both the IBM 370/168 and CDC 6600 machines. A more efficient second version, 
which runs only on the IBM computer, exists. This version contains an 
assembly language subroutine called DINiT with entry points: DREAD, DWRITE, 

DFIND, and DGLOSE. 



C.2 FLOW CHART 


The overall logic for the ARREST computer program is displayed in Figure 
C-1. The main modules of this program comhine all the modal data (Module 1) 
contained on RESIST output tapes onto a single tape and compute the RMS 
deflections and stresses (Module 2). The remainder of the logic is concerned 
with locating the peak tile and primary structure responses for rapid 
identification. 
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SORT 
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jTH TILE c 

AND u*S 






NO. OF EACH TILE 
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Fi^re CM Arrest Program 







.C.3 DATA SET ALLOCATIONS 


The following chart identifies the iinit and file nijunber assignments for 
each data set . Symbols used in the chart, are ; 

6 Deflections 
e Strains 
a Stresses 
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c.ll SUBPROGRAM CALLING SEQUENCE AND COMNiON BLOCKS UTILIZED 

The following table lists the subprograms that are called by each parent 
routine. The conmion blocks used in each routine are also tabularized. The 
subprograms listed under the "Calls” heading appear in the order in which they 
are called. The order of the subprograms listed in the "Subroutine" column 
corresponds to the order in which they appear in the FORTRAN listing. An 
alphabetical index which cross-references this order-of -appearance number is 
provided in Section C.5 



SUBROUTINE 


MAIN 


lODATA 


eOMPAK 


RMS 


SIGDEF 


cabled from 


MAIN 


MAIN 


MAIN 


MAIN 


CALLS 

lODATA 

GETDIM 

GETROW 

DCLOSE 

COMPAK 

RMS 

SIGDEF 

TSTMAX 

BIGBRD 

LBD 

GETDIM 

PUTLAB 

GETROW 

PUTROW 

DCLOSE 

GETDIM 

GETROW 

DCLOSE 

PUTLAB 

PUTROW 

SQRT 

GETDIM 

GETROW 

DCLOSE 

DPRIME 

SORT 

DPMAX 

DTILE 

DTl-JAX 

STILE 

STMAX 


COMMON 

CTAPES 

DIMEN 

FILE 

PROPER 


CTAPES 

FILE 

PROPER 

CTAPES 

FILE 


FILE 

PROPER 


CTAPES 

DIMEN 

FILE 

PROPER 



MO. 

SUBROUTINE 

CALLED FROM 

CALLS 

COMMON 

5 (cont) 



SPRIME 





SPMAX 





STRNGR 


6 

BTILE 

SIGDLF 

- 

CTAPES 

7 

DTMAX 

SIGDEF 

- 

CTAPES 

8 

STILE 

SIGDEF 

- 

CTAPES 

9 

STMAX 

SIGDEF 

- 

CTAPES 

10 

DPRiME 

SIGDEF 

GETDIM 

CTAPES 




GETROW 

DIMEN 




DOLOSE 

FILE 

11 

DPM&X 

SIGDEF 

GETDIM 

CTAPES 




GETROW 

DIMEN 




DCLOSE 

FILE 

12 

SPRIME 

SIGDEF 

GETDIM 

CTAPES 




GETROW 

FILE 




DCLOSE 

DIMEN 

13 

SPMAX 

SIGDEF 

- 

CTAPES 





DIMEN 

ih 

STRNGR 

SIGDEF 


FILE 





DIMEN 



• 


CTAPES 

15 

SORT 

SIGDEF 

- 

DIMEN 

16 

TSTMAX 

MAIN 

- 

CTAPES 





DIMEN 





FILE 

17 

BIGBRD 

lODATA 

- 

CTAPES 

18 

LDB 

I0DATi\ 

EOFOI 

-• 




PLB 


19 

PLB 

LDB 

- 

- 

20 

EOFOI 

LDB 

- 

- 
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NO, 

SUBROUTINE 

GALLED FROM 

CALLS 

COMMON 

21 

PUTLAB 

COMPAK 

DFIND 

CTAPES 



RMS 

DWRITE 

FILE 





PUTGET 

22 

PUTROW 

COMPAK 

DWRITE 

- 



RMS 

PACK 





UNPACK 

• 




DCLOSE 


23 

GETDIM 

MAIN 

DFIND 

ct;^es 



COMPAK 

DREAD 

FILE 



RMS 


PUTGET 



SIGDEF 





DPRIME 





DPMAX 





SPRIME 




GETROW 

MAIN 

DREAD 

- 



COMPAK 

UNPACK 




RMS 

PACK 




SIGDEF 





DPRIME 





DPMAX 




• 

SPRIME 



25 

DWRITE 

PUTLAB 

PUT 

- 



PUTROW 



26 

PUT 

DWRlTE 

- 


27 

BREAD 

GETDIM 

GET 

- 



GETROW 



28 

GET 

DREAD 



29 

DFIND 

PUTLAB 

- 

- 



GETDIM 



30 

DOLOSE 

MAIN 

- 

- 
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NO. 

SUBROUTINE 

CALLED FROM 

CALLS 

COMMON 

30 (cont) 


COMPAK 

RMS 

SIGDEF 

DPRIME 

DPMAX 

SPRIME 



31 

PACK 

PUTROW • 
GETROW 

- 

- 

32 

UNPACK 

PUTROW 

GETROW 

- 

- 
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G,5 BRIEF DESCRIPTION OP SUBROUTINES 

An alphabet ir-al and numerical index, cross-referencing the 32 subroutines 
contained in ARREST, appears in Figure C-2. A brief description of each sub- 
routine is also presented in the order in which it appears in the FORTRAN 
program. 


C-11 




n 

ORDER 

1 

SUBROUTINE 

ORDER 

SUBROUTINE 

ORDER j 

ALPHABETICAL LISTING OF SUBROUTINES 

1 

BIGBRD 

17 

COMPAK 

3 

DOLOSE 

30 

DPI NO 

1 

29 

DPMAX 

11 

DPRIME 

10 

DREAD 

27 

DTILE 

6 

DTMAX 

1 

7 

DWRITE 

25 

EOF01 

20 

GET 

28 

GETDIM 

23 

GETROW 

24 

lODATA 

2 

LOB 

18 

MAIN 

1 

PACK 

31 

PLB 

19 

PUT 

26 

PUTUAB 

21 

PUTROW 

22 

RMS 

4 

SIGDEF 

5 

SORT 

15 

SPMAX 

13 

SPRIME ' 

12 

STILE 

8 

STMAX 

9 

STRNGR 

14 

TSTMAX 

16 



UNPACK 

32 

i 



NUMERICAL INDEX 


MAIN 

1 

IGDATA 

2 

COMPAK 

3 

RMS ^ 

4 

SIGDEF 

5 

DTILE 

6 

DTMAX 

7 

STILE 

8 

STMAX 

9 

DPRIME 

10 

DPMAX 

11 

SPRIME 

12 

SPMAX 

13 

STRNGR 

14 

SORT 

15 

TSTMAX 

16 

BIGBRD 

17 

LDB 

18 

PLB 

19 

EOF01 

20 

PUTLAB 

21 

PUTROW 

22 

GETDIM 

23 

GETROW 

24 

DWRItE 

25 

PUT 

26 

DREAD 1 

27 

GET 

28 

0FINO 

29 

DOLOSE 

30 

PACK 

31 

UNPACK 

32 




Figure C-2 Subroutine Cross-Reference 



C '^12 















DESCRIPTION OP SUBROUTINES 


1 M.'VIN - Initialises input/output unit nimi'bers and calls all input and 
data processing routines. 

2 lODATA - Beads all data from input stream and prints user options. 

3 COMPAK - Incorporates modal data from all tapes onto the tape for the 
first mode in a non-restart job. 

4 IlMS - Generates all RMS deflections, strains and stresses and stores 
them for future printout. 

5 SIGDEF - Supervisor that calls routines for printout of EMS results. 

6 DTILE - Prints RMS TPS displacements for tile. 

T DTMAX - Prints maximum RMS TPS displacements for tile. 

8 STILE - Prints RI-3S stresses for tile. 

9 STfiAX - Prints maxira’jm RMS stresses for tile. 

10 DPRIliE - Prints primary structure EMS deflections. 

11 DPMAX - Prints maximum primary structure EMS deflections. 

12 SPRIIiE - Prints mid-point plate member RI43 strains and stresses. 

13 SP;-IAX: - Prints maximum raid-point plate member RMS strains and stresses. 
lU 3TBNC-E - Prints stringer 3tra.ins and stresses. 

15 SORT - Sorts set of Fi-lS rerults into decreasing order. 

16 TSTL-'!AX ^ Prints maximum EMB streos for each tile. 

17 BIGB3D - Prints AEPEUT title sheet before program output. 

18 LDL - Reads and lists dct-i from input stream and generates data file 

user's program. 

19 PLB - Starts a new page ar.i./or skips a number of lines. 

20 EOFOI - Checks for an end of flic on input data set unit, sets end 
of f -■ le control worn ui.tJon to its proper value and reads a two card 
problem title. 

f-13 
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21 PUTLAB - Puts label on data set containing matrix dimensions and prints 
out label information, if requested. 

22 PUTROW - Writes row of matrix on data set in format designated by packing 
factor . 

23 GETDIM - Gets label of matrix from data set and prints out label 
information, if requested. 

2k GETROW - Reads row of matrix from data set in format designated by 
packing factor. 

25 DWRITE - Converts number of bytes in array into words by dividing by 
four. 

26 PUT - Writes singly-dimensioned array onto data set. 

27 DREAD - Converts number of bytes in array into words by dividing by 
four. 

28 GET - Reads singly-dimensioned array from data set. 

29 DFIND - Rewinds data set to appropriate file. 

30 IKJLOSE - Rewinds data set. 

31 PACK - Packs rows of matrix so they may be written on a data set in an 
efficient manner. This is done by representing strings of zeros by 

a single fixed' point negative integer where the value of the integer 
represents the number of zeroes in the string. Ron-zero nianbers are 
preceded by a fixed point number indicating the number of non-zero 
ntimbers that follow. 

32 UNPACK - Unpacks rows of matrices that have been packed by subroutine 
PACK. 
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C.6 TYPICAL SUBPROGRAM OUTPUT 


The following listings display subprogram output with all the debugging 
clues tiirned on. The name of the particular subprogram which generated each 
set of output is printed in bold type on the rotated right hand side of each 
photographed computer page. 
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